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ABSTRACT;  This  report  describes  the  production  and  use  of  a 
magnetic  tape  library  of  thermodynamic  data.  The  tape  is 
produced  from  a  deck  of  punched  cards,  six  for  each  of  more 
than  one  hundred  chemical  species.  The  data  on  each  species 
Include  an  identification  number,  coded  references,  formula, 
and  atomic  weight  (of  elemental  species),  the  enthalpy  change 
from  O^K  to  29o.l5®K,  the  standard  heat  of  formation  at 
298.15°Kj  melting  point  and  heat  of  fusion,  boiling  point  and 
heat  of  vaporization,  coefficients  giving  the  specific  heat  as 
a  rational  function  of  the  temperature,  and  integration  constants 
permitting  calculation  of  enthalpy  and  entropy.  All  data  are 
given  for  a  pressure  of  one  atmosphere.  These  data  permit  cal¬ 
culation  of  enthalpy,  entropy  and  the  Gibbs  free  energy  over  the 
temperature  range  1000*^  -  6000‘^K,  The  machine  program  used  in 
the  v;ork  is  presented  accompanied  by  illustrations  of  its  use. 
Sources  of  the  original  data  are  given.  Examples  are  produced 
from  the  literature  showing  discrepancies  for  some  aluminum  and 
boron  compounds.  Future  work  will  extend  tape  coverage  to  more 
than  350  species,  allow  for  variation  of  the  form  of  the  speci¬ 
fic  heat  function  and  provide  for-  updating  of  the  data  tape. 


PUBLISHED  SEPTEMBER  1962 


APPROVED  BY;  D  Price,  Acting  Chief 

Physical  Ohemlstry  Division 
CHEMISTRY  RESEARCH  DEPARTMENT 
U.  S.  NAVAL  ORDNANCE  LABORATORY 
WHITE  OAK,  SILVER  SPRING,  MARYLAND 


1 


NOLTR  62-84 


15  May  I962 


This  report  describes  the  construction  of  a  magnetic  tape 
library  of  thermodynamic  data,  for  use  with  a  high  speed  digital 
con5)uter.  The  work  was  carried  out  as  a  part  of  the  Naval 
Ordnance  Laboratory's  propellant  development  studies  under  Task 
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INTRODUCTION 


The  results  presented  In  this  report  are  a  by-product  of 
the  Naval  Ordnance  Laboratory's  effort  in  the  field  of  rocket 
propellant  research.  Although  taped  thermodynamic  data  can  be 
used  for  other^ problems,  the  tape  was  made  specifically  for  use 
with  a  program'"’  developed  here  to  compute  the  theoretical 
performance  parameters  of  solid  rocket  propellants.  As  is  well 
known,  the  computation  of  specific  inpulse  and  other  character¬ 
istic  parameters  of  a  propellant  involves  the  use  of  large 
amounts  of  data  for  many  chemical  species;  to  be  most  useful 
for  machine  computation  these  data  should  be  easily  accessible, 
and  amenable  to  simple  modification  if  Improved  values  are 
received.  With  the  present  state  of  the  computer  art,  the  use 
of  a  magnetic  tape  best  fulfills  the  above  requirements. 
Accordingly  the  objective  was  to  collect,  verify  and  collate 
thermodynamic  data,  and  store  them  on  a  magnetic  tape  in  a 
practical,  accessible  form. 

For  making  the  tape,  the  data  fall  into  two  categories: 
part,  such  as  heats  of  formation,  can  be  read  onto  the  tape 
immediately  after  verification,  while  the  remainder,  such  as 
specific  heats,  must  be  reduced  to  a  form  suitable  for  tape 
storage.  Early  workers  tended  to  store  data  tables,  giving 
values  at  discrete  temperature  values.  Such  data  required 
Interpolation  for  Intermediate  temperatures.  A  better  practice 
is  to  reduce  the  data  by  a  curve  fitting  process  to  a  rational 
function  of  the  temperature  and  store  only  the  coefficients  of 
the  function.  A  Fortran  language  program,  to  be  discussed 
later,  was  written  and  compiled  to  do  this  as  well  as  other 
parts  of  the  data  work.  After  collection  and  preparation,  the 
data  were  punched  on  cards  and  subsequently  written  onto  an  lEM 
magnetic  tape.  In  this  form  they  may  be  readily  recalled  for 
further  work  by  the  IBM  704  or  7090  computers.  Sample  Data 
Program  output  sheets  are  illustrated  £ind  discussed  in  the  body 
of  the  report. 

The  various  data  sources  searched  agree  with  each  other 
within  one  percent  for  the  common  compounds  of  carbon,  hydrogen, 
nitrogen  and  oxygen  and  the  common  halides.  In  some  cases, 
particularly  with  newer  types  of  propellant  ingredients  e.g.. 


*This  program  has  been  in  operation  since  March  i960.  It  will 
be  described  in  a  forthcoming  report. 
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light  metals,  the  necessary  data  are  not  in  final  form  as 
Bhovm  by  the  lack  of  agreement  among  competent  workers  in  the 
field.  Again,  evidence  of  new  species  stable  at  high  tempera¬ 
tures  is  being  found  and  for  these  the  data  are  still  in  the 
preliminary  stage.  More  refined  studies  must  be  made  on  some 
ternary  and  higher  compounds. 

Although  a  data  sheet  has  been  printed  out  for  verifica¬ 
tion  for  each  species  on  the  tape,  it  is  not  the  intention  here 
to  add  another  member  to  the  already  large,  and  still  growing 
list  of  published  collections  of  thermodynamic  data.  Since 
most  of  these  lists  go  back  to  a  relatively  few  common  sources 
such  as  reference  (l),  they  represent  duplication  of  effort  in 
the  case  of  the  well  known  species  and  multiplication  of 
inaccurracles  in  the  case  of  newer  ones.  The  latter  require 
further  effort  to  eliminate  Inaccurracles;  discussion  and 
selection  of  such  data  are  better  Included  in  reports  of  work 
utilizing  the  particular  quantities.  For  example,  reported 
theoretical  specific  impulses  should  always  be  accompanied  by 
input  data  on  any  questionable  species  considered  in  the  calcu¬ 
lations  . 

In  the  meantime,  since  there  are  differences  between  one 
laboratory  and  another  in  the  calculatlonal  methods  used,  it 
would  be  helpful  if  all  groups  at  least  used  common  data.  The 
JANAF  Thermo chemical  Panel  has  undertaken  the  ambitious  task  of 

bringing  this  about;  and  one  of  the  final  objectives  of  the 

present  work  is  to  bring  the  data  library  tape  into  accord  with 
the  JANAF  Thermochemlcal  Tables.  Further,  coverage  is  being 
extended  from  105  to  about  550  chemical  species  and  the  program 

is  being  modified  to  allow  more  latitude  in  the  treatment  of 

specific  heats  and  other  data. 

THEORy 

As  noted  In  Section  I,  the  data  to  be  recorded  were  divided 
into  two  parts;  those  to  be  taped  Immediately  and  those  to  be 
reduced  first  to  a  set  of  storable  parameters.  They  were  assem¬ 
bled  mainly  from  the  sources  given  by  references  1  to  11.  Data 
such  as  heats  of  formation,  melting  and  boiling  points  were 
collected,  verified,  punched  on  IBM  cards  and  read  directly  onto 
the  magnetic  tape. 

For  data  from  which  thermodynamic  functions  were  to  be 
derived,  the  procedure  was  different;  the  items  to  be  stored  on 


2 


NOLTR  62-84 


tape  for  each  species  were:  coefficients  of  a  function  of  the 
temperature,  T,  which  defined  the  specific  heat,  C  ,  at  a  , 
pressure  of  one  atmosphere,  and  two  Integration  constants,  H 
and  S',  discussed  below.  Prom  these  data  by  the  application  of 
elementary  thermodynamics,  the  enthalpy  H,  the  entropy  S  and 
the  Gibbs  free  energy  0  were  derived. 

Since  Cp,  an  experimental  quantity  put  into  an  appropriate 
analytical  fdrra,  has  to  be  integrated  to  obtain  the  enthalpy,  it 
has  to  be  represented  by  some  function  of  Tj  and  no  one  function 
provides  a  representation  equally  good  for  all  chemical  species. 
As  pointed  out  by  Maler  and  Kelley  (12)  it  was  at  one  time  com¬ 
mon  practice  to  represent  the  enthalpy  by  an  elementary  power 
series  in  T,  terms  as  high  as  usually  being  taken.  Differen¬ 
tiation  then  resulted  in  a  quadratic  function  for  Cp  which 
could,  " . . . .  have  a  maximum  in  the  temperature  range  under  con¬ 
sideration  due  to  the  signs  and  magnitudes  of  the  coefficients, 
but  such  a  maximum  does  not  correspond  to  known  physical  facts". 
They  proposed  the  addition  of  negative  powers  of  the  temperature, 
and  showed  by  comparative  results  the  ln5)rovement  obtained  over 
the  simple  quadratic  function.  Since  then  several  different 
functions  have  been  eriployed,  with  varying  success,  to  represent 
Cp(T)  analytically.  The  form  used  here  was  noted  in  a  Wright 
Air  Development  Center  report  (15),  and  has  since  been  used  in 
the  preparation  of  the  JANAF  Thermocheralcal  Tables  (21),  For 
all  species  they  use  the  function 

Cp  =  Cp(T)  =  a  +  be  +  c-e-  ^  +  d  -e-  ^  +  e^  (l) 

where,  for  convenience 

■©•=  10"^  T  (2) 

The  problem  was  then  to  determine  the  set  of  parameters  a, 
b,  c,  d  and  e  for  each  chemical  species.  Obviously  five  equa¬ 
tions  were  needed  and  these  were  obtained  by  using  experimental 
values  of  Cp  at  five  different  tenperatures.  There  resulted  a 
set  of  flve^equatlons  linear  in  a,  b,  c,  d  and  e,  and  these 
were  solved  readily  using  a  single  matrix  inversion  technique. 

With  a,  b,  c,  d  and  e  known  for  a  species,  Cp  was  defined; 
for  the  enthalpy  and  entropy,  it  remained  only  to^determine  the 
integration  constants  H'  and  S'  respectively.  The  elementary 
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steps  needed  are  summarized  in  Appendix  and  a  detailed 
example  was  presented  in  reference  (15).  The  variation  of  H' 
and  S'  was  shown  to  be  approximately  1^  between  1000°  and 
6000°K.  This  was  considered  tolerable. 


DATA  COLLECTION 

Although  several  score  sources  have  been  examined  in  the 
process  of  collecting  thermodynamic  data,  only  a  few  of  the  more 
Inportant  ones  are  referenced  here  (1  -  11).  Of  these,  the  most 
In^ortant  were  from  the  National  Bureau  of  Standards  (l,  4,  6, 
7).  The  reference  states  of  the  elements  were  taken  as  given  by 
Stull  and  Sinke  (lO). 

Among  the  data  in  use,  good  agreement  exists  between  one 
laboratory  and  another  for  the  common  compounds  of  carbon, 
hydrogen,  nitrogen  and  oxygen,  and  also  of  fluorine  and 
chlorine.  However,  when  metallic  combustion  products  are  con¬ 
sidered,  agreement  is  often  lacking  between  one  laboratory  and 
another.  Present  knowledge  of  these  compounds  has  been  summar¬ 
ized  by  Dobbins  (9),  pages  99  et  seq;  and  Ollles  (l4)  has 
reviewed  the  status  of  work  on  high  temperature  ternary  species, 
another  area  where  general  agreement  is  lacking.  It  is  necess¬ 
ary  here  only  to  present  a  few  illustrative  examples  of  varia¬ 
tion  in  data  used. 

•ic  . 

Pig.  1  shows  three  derived  curves  of  common  log  K  vs  1/T 
for  the  reaction 

AlgO^Cl)  +  2H2(g)<_->  Al20(g)  +  2H20(l). 

The  upper  and  lower  curves  represented  the  extreme  values  among 
those  based  upon  material  in  several  classified  and  unclassified 
sources  (2,  l6,  17,  l8,  19,  20).  Since  the  scale  on  the  ordin¬ 
ate  represented  powers  of  ten,  the  graph  Illustrated  the  large 
differences  from  one  laboratory  to  another.  The  third  curve 
shown  was  based  upon  values  vrtilch  are  now  considered  acceptable, 
and  are  reported  in  the  open  literature  (21). 


^Before  the  publication  of  reference  (15)j  flame  temperatures 
and  compositions  were  computed  here  using  the  older  equili¬ 
brium  constant  method,  thus  the  Interest  in  log  K.  Although 
the  method  is  no  longer  used,  the  results  shown  in  Pig.  1  are 
typical  of  the  kind  of  disagreement  often  found. 
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Another  exarrple  Is  the  standard  heat  of  formation  of  boron 
nitride.  This  has  been  variously  reported  (in  kcals  per  mol)  as 
131.4  (5),  155  (9)  and  158. 1  (l).  Again  the  existence  of  a 
boron  subchloride  BCl^  has  been  inferred  (9)  from  the  existence 
of  the  dimer;  and  the  neat  of  formation  of  the  subchloride  has 
been  reported  as  -24.6  (5)  and  -82.0  (9). 

In  cases  such  as  these,  the  consistency  of  the  given  data 
with  other  known  values  for  the  species  was  checked.  Until 
recently  (See  Recommendations)  values  given  by  the  National 
Bureau  of  Standards  were  considered  to  be  the  most  acceptable  in 
cases  which  could  not  otherwise  be  decided.  Henceforth,  values 
from,  the  JANAP  tables  will  be  used. 


DATA  CARD  DECK 

Six  IBM  cards  are  used  to  hold  the  data  for  each  species 
prior  to  taping  them.  The  first  three  cards  (of  the  set  of  six) 
contain  thermodynamic  data  and  are  described  In  Pigs,  2,  3  and 
4.  Fig.  5  typifies  the  last  three  cards  of  the  set.  The 
summary  of  data  for  one  species  is  indicated  in  Table  1.  In 
this  table  the  lower  case  letters,  '’cc"  in  the  left  hand  column 
stand  for,  "card  column"  and  define  the  card  fields.  The  center 
column  of  the  table  lists  the  Information  in  the  corresponding 
card  columns.  The  column  on  the  right  side  of  the  table  gives 
the  storage  designations  for  the  data  when  they  are  being  used. 
The  letters  SPD  in  this  column  originally  meant  Species 
Data,  and  are  now  used  to  identify  all  data  storage  designations. 

The  data  for  the  first  three  cards  of  a  set  may  be  punched 
as  soon  as  verified.  The  last  three  cards,  typified  by  Pig,  5, 
refer  to  the  gas  (card  4),  liquid  (card  5)  and  solid  (card  6) 
phases  of  a  species.  Each  of  these  three  contains  seven 
numbers.  The  first  five  are  respectively  the  coefficients  a,  b, 
c,  d  and  e  for  the  species  and  the  last  two  are  the  correspond¬ 
ing  Integration  constants  H’  and  S',  for  the  enthalpy  and 
entropy  respectively. 


MACHINE  PROGRAM 

On  page  three  and  Appendix  A  equations  are  presented  for 
the  specific  heat  at  constant  pressure,  the  enthalpy,  the  entro¬ 
py  and  Gibbs  free  energy  as  functions  of  temperature.  A  Fortran 
program  was  written  to  derive  the  coefficients  of  the  Cp  equa¬ 
tion  and  to  compute  tables  of  these  functions.  A  block  diagram 


5 


NOLTR  62-8^+ 


Of  the  program  Is  shown  In  Pig.  6  and  a  copy  of  the  Fortran 
listing  is  Included  as  Appendix  B. 

Five  options  were  built  into  the  program.  Including  one 
"dummy"  for  future  use.  The  options  have  extended  the  program's 
usefulness.  According  to  the  Job  to  be  done,  the  proper  option 
number  was  punched  on  a  card  and  the  program  followed  the  indi¬ 
cated  path.  A  description  of  the  various  options  follows,  in 
logical  rather  than  numerical  order. 

Option  (2).  When  new  experimental  Cp  data  were  acquired, 
we  wished  to  find  the  coefficients  of  the  Cp  function.  The 
cards  were  punched  as  shown  in  Pig.  7.  This  option  caused  the 
program  to  coir5)ute  a,  b,  c,  d,  e,  H'  and  S',  and  print-out  the 
results  for  inspection..  Figure  8  is  a  sample  of  the  print-out 
for  this  option  showing  the  input  Cp,  H  and  S  and  the  resulting 
set  of  constcints. 

Option  (-1).  If  the  data  table  derived  from  the  above 
computed  coefficients  was  desired,  the  cards  for  this  option 
were  punched  as  shown  in  Pig.  9  and  used  in  conjunction  with 
option  (2).  A  sample  of  this  output  is  shown  in  Pig.  10.  The 
machine  time  for  the  exan5)le  shown  was  approximately  five 
seconds  on  the  IBM  704. 

Now  if  examination  showed  an  error  in  the  above  table 
(Pig.  10)  it  was  corrected  at  this  point.  If,  however,  the 
data  were  consistent  with  other  known  data,  the  new  coefficients 
were  punched  on  card  No.  4,  5  or  6  for  this  species,  depending 
upon  the  phase  being  considered  and  were  now  ready  to  be  read 
onto  the  tape.  In  this  way,  the  tape  library  was  compiled. 

The  tape  is  kept  in  the  computer  room,  and  the  corresponding 
card  deck  Is  kept  for  use  in  future  modifications  of  the  tape. 

Figures  11  and  12  illustrate  the  arrangement  of  the  data 
on  tape.  It  is  exactly  as  the  data  card  deck  is  read  in. 

Fig.  11  shows  the  first  cards  (of  each  set  of  six)  for  several 
species.  The  identification  numbers  are  in  the  left  hand 
column,  and  the  first  forty  rows  were  reserved  for  atomic 
species.  A  figure  one  occurs  in  the  column  corresponding  to 
the  identification  number  of  the  atomic  species,  l.e.,  row  12 
has  a  one  printed  in  column  12,  row  l4,  a  one  printed  in  column 
l4,  etc.  This  permitted  a  simple  method  of  defining  molecular 
types,  which  occupy  the  rows  beyond  number  forty.  For  example, 
a  one  in  row  7  identifies  atomic  hydrogen  and  a  one  in  row  11 
identifies  atomic  oxygen.  So  in  order  to  specify  a  molecule  of 
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HgO  (row  45),  a  2  was  placed  in  column  7  and  a  1  in  column  11. 

The  assignment  of  a  column  to  a  particular  element  allowed 
efficient  formulation  of  molecular  species  on  tape.  Upon  recall 
of  the  data  a  Symbol  Sub-routine  identified  the  numbers  in  each 
row  and  converted  them  to  atomic  formulae.  The  remaining  infor¬ 
mation  (right  hand  columns  of  Pig.  12),  includes  the  weights, 
given  for  atomic  species  only,  and  a  set  of  numbers  giving  the 
references  for  each  species. 

Pig.  12  is  a  section  of  the  print-out  corresponding  to  the 
last  five  data  cards  for  each  of  nine  separate  species.  What 
appears  are  three  rows  for  each  species,  corresponding  to  data 
cards  2,  3  and  4.  This  is  because  the  species  here  listed  do 
not  have  liquid  (card  5)  or  solid  (card  6)  phases  at  the  temper¬ 
atures  being  considered.  The  full  row  in  each  group,  shows  in 
order  a,  b,  c,  d,  e,  H'  and  S*  for  the  gas  phase.  If  liquid  or 
solid  phases  for  any  of  these  species  were  present,  they  would 
be  placed  in  the  two  rows  immediately  below  that  for  the  gaseous 
phase.  This  option  has  been  useful  when  a  problem  result  indi¬ 
cated  a  possible  error  on  the  tape.  It  provides  for  fast 
checking , 

Option  (1).  When  it  was  desired  to  examine  the  tape 
contents,  the  input  cards  were  punched  as  in  Pig.  I3.  The 
print-out  from  this  option  consists  of  the  taped  information 
discussed  above,  listed  with  proper  English  column  headings. 

It  has  been  a  convenient  means  of  examining  the  tape  contents. 

Option  (3).  As  in  the  other  options,  the  cards  used  here 
(Pig.  l4)  were  slnple  to  prepare.  The  data  tables  for  any 
species  (Plg.  10),  or  chosen  group  of  species,  were  obtained  as 
the  specific  Inpulse  program  uses  them.  When  desired,  this 
option  can  print  data  tables  for  all  the  species  on  tape;  such 
tables  are  in  suitable  form  for  publication. 

In  using  the  program,  data  cards  are  punched  according  to 
the  appropriate  figures  for  each  option.  The  cards  are  then 
placed  at  the  end  of  the  program  deck  and  run  into  the  card 
hopper.  Obviously,  Options  T-l)  and  (2)  operate  directly  from 
card  input,  while  (l)  and  (3)  require  that  a  data  tape  reel  be 
placed  on  tape  unit  No.  7,  the  data  input  unit  of  the  704  com¬ 
puter,  For  the  IBM  709O,  the  tape  number  is  A7,  If  by  the  use 
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of  Options  (-1)  and  (2),  acceptable  species  data  have  been 
generated,  they  are  punched  onto  card  4,  5  or  6,  and  Inserted 
Into  the  library  deck  along  with  the  first  three  cards  for  the 
species.  In  this  way  data  on  any  species  In  the  library  deck 
can  be  replaced. 

The  library  deck  and  the  tape  are  the  end  results  of  the 
work.  The  program  deck,  used  to  handle  the  data,  may  be  repro¬ 
duced  simply  by  punching  the  cards  shown  In  the  Fortran  listing 
of  Appendix  B.  After  this  program  deck  Is  compiled.  It  Is 
necessary  only  to  choose  an  option  In  order  to  use  the  system, 
Table  2  Is  a  list  of  the  103  species  for  which  data  are  on  tape. 

RESULTS 

A  continuing  system  has  been  set  up  for  the  acquisition, 
verification  and  recording  on  magnetic  tape,  of  thermodynamic 
data.  Methods  have  been  developed  for  replacing  data  now  on 
tape,  and  for  adding  data  on  new  species.  This  work  has  been 
carried  on  simultaneously  with  work  on  the  theoretical  pro¬ 
pellant  calculation  program  so  that  a  modicum  of  data  was 
available  as  soon  as  the  latter  was  ready  to  be  tested.  Much 
Information  has  since  been  added  to  the  library,  and  It  has  been 
used  successfully  for  hundreds  of  propellant  calculations.  The 
Laboratory  has  not  essayed  the  publication  of  a  data  book 
because  the  already  long  list  of  such  books  precludes  further 
duplication.  However,  a  prlnt-out  of  the  material  on  the  tape 
has  been  made  for  our  own  use,  and  the  latest  tape  revision 
occurred  In  May  196I. 

It  may  be  noted  that  the  present  method  of  storing  thermo¬ 
dynamic  data  Is  a  considerable  improvement  over  older  methods. 
Stored  data  tables  are  too  unwieldy  since  it  Is  at  times  desir¬ 
able  in  a  calculation  to  consider  as  many  as  20  -  30  chemical 
species  and  the  necessary  table  storage  space  Is  too  large, 
even  for  modem  memory  devices.  Besides  this,  the  storage  of 
data  tables  requires  the  use  of  Interpolation  subroutines  since 
the  occurrence  of  a  temperature  value  lying  between  two  of 
those  given  In  the  table  Is  inevitable. 


FUTURE  WORK  AND  RECOMMENDATIONS 

A  better  method  of  curve  fitting  might  be  used  in  the 
Fortran  program.  A  least  squares  fit  and  Tchebicheff  Polynom¬ 
ials  are  both  being  considered.  Certain  techniques  which  here¬ 
tofore  have  been  separate  from  the  data  program  should  be 
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assimilated  into  it  as  further  options.  For  example,  one 
separate  program  computes  the  coefficients  a,  b,  c,  d  and  e  in 
equation  (l)  from  coefficients  as  given  by  other  Cp  equations 
for  the  same  species;  while  another  program  prints  out  all  the 
data  available  at  any  given  temperature,  as  a  single  block  of 
Information,  Further,  techniques  are  now  available  for  provid¬ 
ing  relatively  automatic  updating  of  the  tape.  These  techniques 
would  allow  an  even  sin^ler  method  of  tape  modification. 

Since  many  species  of  great  Importance  in  rocket  propellant 
calculations  are  still  not  completely  defined,  a  continuing  data 
search  should  be  made  for  the  most  probable  values  among  those 
available  for  these  species. 

In  the  past,  the  National  Bureau  of  Standards  data  has 
been  used  as  the  most  reliable.  In  the  interest  of  furthering 
consistency  among  all  groups  in  the  field,  the  data  used  by 
this  Laboratory  should  be  based  upon  the  data  published  by  the 
JANAF  Thermochemical  Panel  (21,  22).  Revisions  of  these  data 
are  continually  being  supplied  by  the  Panel;  so  the  fact  that 
the  data  in  many  cases  is  admittedly  tentative  should  not  be  a 
deterrent . 

Finally  it  would  be  of  considerable  Interest  to  extend  this 
work  in  the  direction  of  more  extreme  thermodynamic  conditions 
because  of  its  potential  use  in  explosives  confutations. 
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TABLE  1 

Arrangement  of  Species  Data 


Card 

Card  Columns 

Data 

Tape 

I 

cc 

1 

...  3 

ID  Number 

SPD  (1,1) 

cc 

4 

..  43 

Formulation 

Not 

Stored 

cc 

44 

..  53  , 

Mol.  Wt. 

SPD 

(1,2) 

cc 

54 

{0.1) 

Phase 

SPD 

1,3 

cc 

55 

(I,-.) 

Ref. 

Not  Stored 

II 

cc 

1 

..  10 

(g) 

SPD(I,4) 

cc 

11 

..  20 

(1) 

SPD(I,5) 

cc 

21 

..  30 

( B ) 

spd(i,6) 

cc 

31 

..  40 

bp 

SPD(I,7) 

cc 

41 

..  50 

top 

SPD(I,8) 

III 

cc 

1 

..  10 

H2^8“^o  {^^ 

SPD(I,9) 

cc 

11 

..  20 

H298-H0  (1) 

SPD(I,10) 

cc 

21 

..  30 

H298-H0  (s) 

SPD(1,11) 

cc 

31 

..  40 

AHv 

SPD(I,12) 

(vaporization) 

cc 

41 

..  50 

AHf( fusion) 

SPD(I,13) 

IV 

cc 

1 

..  10 

a(g) 

SPDI 

1,14 

cc 

11 

..  20 

b 

SPDI 

1.15 

cc 

21 

..  30 

c 

SPDI 

I,  lb. 

cc 

31 

..  40 

d 

SPDI 

1,17 

cc 

41 

..  50 

ei 

SPDI 

1,18 

cc 

51 

..  60 

IH' 

SPDI 

1,19 

cc 

6l 

..  70 

2S' 

SPDI 

(l,20| 

V 

cc 

1 

..  10 

a(l) 

SPDI 

[1,21] 

cc 

11 

..  20 

b 

SPDI 

1,22 

cc 

21 

..  30 

c 

SPDI 

1,23 

cc 

31 

..  40 

d 

SPDI 

1,24 

cc 

4l 

..  50 

e 

SPDI 

1,25 

cc 

51 

..  60 

IH' 

SPDI 

1,26 

cc 

6l 

70 

2S' 

SPDI 

[1,27] 

VI 

cc 

1 

..  10 

a(s) 

SPDI 

[1,28] 

cc 

11 

..  20 

b 

SPDI 

1,29 

cc 

21 

..  30 

c 

SPDI 

1,30 

cc 

31 

..  4o 

d 

SPDI 

1,31 

cc 

4l 

..  50 

e 

SPDI 

1,32 

1 

cc 

51 

. .  60 

H' 

SPDI 

1,33 

1 

cc 

6l 

..  70 

S' 

SPDI 

ll,34i 

1 
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TABLE  2 

List  of  Species  on  Tape 


A1 

HF 

CP 

BCl 

LiCl 

Al(l) 

CI2 

CFg 

BClg 

LiP 

B 

GIF 

CF3 

BCl^ 

LiH 

Be 

HCl 

CP4 

BF 

LIO 

C 

CIO 

C2P2 

BF 

LigO 

0(s) 

^2 

CF^H 

BH 

LiOH 

01 

HS 

CHgPg 

BN 

AlgO 

F 

SO 

CH^F 

BN(s) 

AIF^ 

H 

N2 

COFg 

BO 

AlPg 

Li 

NFO 

CH 

AlFO 

Mg 

NH 

CHg 

BgOj 

AlCl^ 

N 

NO 

CH4 

AlP 

AlCl^(l) 

0 

^2 

CN 

AlH 

AlClg 

P 

^4 

C2N2 

AlO 

AICIO 

S 

PH 

CNF 

Al20j(l) 

AlClPg 

Si 

PH^ 

CO 

Al20j(s) 

AICIF 

®2 

PN 

COg 

BeP  ® 

AlCl 

PO 

COHp 

BeH 

AlC 

H2 

^2 

CP 

BeO 

AlP 

OH 

^3 

cs 

MgP 

AIS 

H2O 

CCl 

SIO 

MgH 

^2 

COClg 

®2 

CVI 

•H 

NOTE:  All  species  are  gaseous  unless  noted,  in  which  cases  an 
s  or  1  in  parenthesis  Indicates  the  solid  or  liquid  phase 
respectively. 
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log  K  VS  '/jFOR  THE  REACTION  OF  Al203(^)AND  H2(^) 


FIG.  3  INFORMATION  ON  DATA  CARD  NO. 
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FIG.  4  INFORMATION  ON  DATA  CARD  NO. 
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FIG.  5  INFORMATION  ON  DATA  CARDS  NOS.  4,  5  AND  6 


NOTE!  CURRENT  TAPE  NUMBER  WILL  VARY 


NOLTR  62-84 


THREE  CARDS  ARE  USED  FOR  THIS  OPTION 


Card  No.  1: 

cc  2,  punch:  2 
Card  No.  2: 

punch:  coefficient  n  determination  a  for  H2O  ^  formula 

(atomic  formula  of  species) 


Card  No.  3: 

cc  1-10 

punch : 

Cp{1000) 

cc  11-20 

punch : 

Cp(l600) 

cc  21-30 

punch : 

Cp(3000) 

cc  31-40 

punch: 

Cp(450O) 

cc  41-50 

punch : 

Cp(5700) 

cc  51-60 

punch : 

H'(JOOO) 

cc  61-70 

punch : 

S'(JOOO) 

Note:  On  Card  No.  2  A  represents  a  blank,  cc  means  card 
column. 

Fixed  point;  decimal  point  must  be  punched. 


Figure  7  -  Input  Cards  for  Option  (2) 
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CO 
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FIG. 8  OUTPUT  FROM  OPTION 
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TWO  CARDS  ARE  USED  FOR  THIS  OPTION 


Card  No.  1 


cc 

1,  2 

punch;  -1 

Card 

No.  2 

cc 

1-10 

punch; 

^1 

cc 

11-20 

punch: 

cc 

21-30 

punch ; 

aT 

NOTE;  T^  and  T^  are  the  lower  and  upper  limits  of  the 

temperature  interval  over  which  data  are  desired. 
“T  is  the  size  of  the  step  from  T^  to  T^^  ^ 


NOTE;  Fixed  point;  decimal  point  must  be  punched. 


Figure  9  -  Input  Cards  for  Option  (-1) 
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THERMODYNAMIC  DATA  FOR  H^O  GAS 


CP 

(H-Ho)/T 

S 

(G-Go)/T 

T 

(CALS/ 

(CALS/ 

(CALS/ 

(CALS/ 

> 

o 

n)  MOL  °K) 

MOL  °K) 

MOL  1 

0 

MOL°K 

) 

1000. 

9.85026E  00 

8.57754E  00 

5.55906E 

01 

-4.70131E 

01 

1100. 

1.01561E  01 

8. 70725 E  00 

5.65149E 

01 

-4.78077E 

01 

1200. 

1.04474E  01 

8.84023E  00 

5.73839E 

01 

-4.85437E 

01 

1300. 

1.07229E  01 

8.97455E  00 

5.82045E 

01 

-4.92299E 

01 

1400. 

1.09823E  01 

9.10880E  00 

5. 89822 E 

01 

-4.98734E 

01 

1500. 

1.12256E  01 

9.24190E  00 

5.97214E 

01 

-5. 04795 E 

01 

1600. 

1.14530E  01 

9.37307E  00 

6.04258E 

01 

-5.10527E 

01 

1700. 

1.16649E  01 

9.50172E  00 

6.10983E 

01 

-5.15966E 

01 

1800. 

1.18619E  01 

9.62744E  00 

6.17414E 

01 

-5.21 140E 

01 

1900. 

1.20444E  01 

9. 74991 E  00 

6.23574E 

01 

-5.26075E 

01 

2000. 

1.22132E  01 

9. 86891 E  00 

6. 29481 E 

01 

-5.30792E 

01 

2100. 

1.23689E  01 

9.98430E  00 

6.35151E 

01 

-5.35308E 

01 

2200. 

1.25122E  01 

1.00960E  01 

6.40598E 

01 

-5.39639E 

01 

2300. 

1.26439E  01 

1.02039E  01 

6.45837E 

01 

-5.43797E 

01 

2400. 

1.27646E  01 

1.03082E  01 

6.50877E 

01 

-5,47796E 

01 

2500. 

1.28752E  01 

1.04087E  01 

6.55731E 

01 

-5.51644E 

01 

2600. 

1.29765E  01 

1.05055E  01 

6<60406E 

01 

-5. 5535 IE 

01 

2700. 

1. 30691 E  01 

1.05988E  01 

6.64913E 

01 

-5.58926E 

01 

2800. 

1.31539E  01 

1.06885E  01 

6.69259E 

01 

-5.62374E 

01 

2900. 

1.32316E  01 

1.07749E  01 

6.73452E 

01 

-5.65703E 

01 

3000. 

1.33032E  lOl 

1.08580E  01 

6.77499E 

01 

-5.68919E 

01 

3100. 

1.33692E  01 

1.09380E  01 

6.81406E 

01 

-5.72026E 

01 

3200. 

1.34306E  01 

1.10149E  01 

6.85180E 

01 

-5. 7503 IE 

01 

3300. 

1. 34881 E  01 

1.10890E  01 

6.88826E 

01 

-5.77937E 

01 

3400. 

1.35425E  01 

1.1I604E  01 

6.92351E 

01 

-5.80748E 

01 

3500. 

1.35947E  01 

1.12292E  01 

6.95760E 

01 

-5.83468E 

01 

3600. 

1.36453E  01 

1.12956E  01 

6.99058E 

01 

-5.86102E 

01 

3700. 

1.36954E  01 

1.13598E  01 

7.02249E 

01 

-5. 8865  IE 

01 

3800. 

1.37455E  01 

1.14219E  01 

7.05339B 

01 

-5.91 120E 

01 

3900. 

1.37966E  01 

1.14821E  01 

7.08332E 

01 

-5.93511E 

01 

4000. 

1.38494E  01 

1.15406E  Cl 

7.11233E 

01 

-5.95826E 

01 

4100. 

1.39047E  01 

1.15976E  01 

7. 14045E 

01 

-5.98069E 

01 

4200. 

1.39634E  01 

1.16532E  01 

7.16774E 

01 

-6.00242E 

01 

4300. 

1.40263E  01 

1,17077E  01 

7.19424E 

01 

-6.02347E 

01 

4400. 

1. 40941 E  01 

1.17611E  01 

7.21997E 

01 

-6.04386E 

01 

4500. 

1,41676E  01 

1.18138E  01 

7.24499E 

01 

-6. 06361 E 

01 

NOTE: 

Given  the  necessary  coefficients, 

compute  and 

print 

data  tab 

1  e 

0 

for  species  nomed,  from  1  000  K  to 

6000°K. 

FIG.  10 

PARTIAL  OUTPUT  FROM  OPTION 

(- 

•1) 
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SPECIE 
NUMBER 
0011 
002  1 
003  1 

004  1 

005 
006 
007 
008 
009 
010 
on 
012 
013 
014 
015 
016 
017 
018 
019 
020 
021 
022 
023 
024 
025 
026 
027 
028 
029 
030 
031 
032 
033 
034 
0351 
036 
037 
038 
039 
040 
041 
042 
043 
044 
045 
046 
047 
048 
049 
050 
051 
052 


ATOMIC 

REFS. 

WEIGHTS 

26.980000  8016  Al 

10.8200 

8143  B 

9.0130 

8131  Be 

12.010000  7082  C 

35.4570 

8020  Cl 

19.000000  8210  F 

■1.0080000  8211  H 

6.9400 

8118  Li 

24.32 

7082  Mg 

14.0080 

8208  N 

■16.0000 

8213  0 

30.9750 

7067  P 

32.066 

8003  S 

28.0000 

8002  Si 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 


2 

3 

2 

1  1 
2f  V 

2 
1  I 

2 
11 
1  1 

1  1 

2 


1 

1  • 

1 

1 

1 

1 

1  . 

1 

1 

1 

1 

1 

1 

1 


8018 

8212 

8018 

8018 

7082005 

8019 

7082094 

8016 


FIG.  II  INPUT  TAPE  LISTING  (A) 
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CARD 

NO. 

2  121600, 

3  2250. 


4 

8.91979 

0.04404 

-0.00241 

0.00017 

-0.19660 

-0.9700 

63.8156 

5 

6 

2 

-53300. 

264.7 

3 

3283. 

4 

18.5699 

1.01499 

-.28713 

.02781 

-.73558 

-4.3637 

88.2420 

5 

6 

2 

74500. 

2 

2170. 

4 

8.59423 

0.36810 

-0.08798 

0.00895 

-0.34429 

-1.3240 

59.8720 

5 

6 

2 

-30000. 

3 

2467. 

4 

13.52244 

0.24640 

-0.05383 

0.00372 

-0.87872 

-3.9585 

70.2310 

5 

6 

O 

-120500. 

3 

2794. 

4 

18.97333 

0.70465 

-0.20308 

0.02012 

-1.28501 

-6,4394 

79.8664 

5 

6 

2 

-218000. 

3 

3045, 

4 

24.72236 

0.87981 

-0.25485 

0.02539 

-1.74270 

-9.2157 

84.9078 

5 

5 

2 

-51300. 

3 

3293. 

4 

19.13848 

2.55046 

-0.37588 

0.03702 

-1.04807 

-5.2626 

78.9155 

5 

6 

2 

-162600. 

3 

2765. 

4 

21,26435 

3.33827 

-0.91481 

0.08746 

-2.04527 

-0.1818 

78.9034 

5 

6 

2 

-1.05500, 

3 

2556. 

4 

17.92884 

5.77349 

-1.58420 

0.15200  ■ 

-2.30013 

-8.8788 

70.8631 

5 

6 

NOTE:  Rows  numbered  2  contain  A  HT(g)  for  each  species. 

Rows  numbered  3  contain  (H290“Ho)  (g)  for  each  species. 
Rows  numbered  4  contain  a,D,c,d,e,  H'and  S'  for  gas  phase. 
Rows  numbered  5  and  6  are  blank,  since  cniy  gas  phases  are 
important  for  the  above  species. 

FIG.  12  INPUT  TAPE  LISTING  (B) 
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TWO  CARDS  ARE  USED  FOR  THIS  OPTION 

Card  No.  1 
cc  2  punch:  1 

Card  No.  2 

Start  in  cc  1,  and  punch  as  follows: 

XX  U  XX  a  XX 

month  Day  Year 

NOTE:  a  on  Card  No.  2  represents  a  blank 


Figure  15  -  Input  for  Option  (1) 
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TWO  CARDS  ARE  USED  FOR  THIS  OPTION 

Card  No.  1 
cc  2  punch:  3 
Card  No.  2 


cc  1,  2,  3 

punch: 

Identification  number 

cc  4 

punch: 

g,  1,  or  8  (phase) 

cc  5,  14 

punch: 

*1 

cc  15,  24 

punch: 

cc  25,  34 

punch : 

AT 

NOTE:  Fixed  point:  decimal  point  must  be  punched. 


Figure  14  -  Input  Cards  for  Option  (3) 
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APPENDIX  A 

Equations  Used  by  the  Program 


Equation  (l),  Section  II,  defines  Cp(T)  at  a  pressure  of 
one  atmosphere;  and  from  Equation  (2) 


dT  =  lO^d^  . 

The  enthalpy  Is  then  obtained  by  Integrating 

dH  =  C  dT  =  10 Vd^  ;  (Al) 

thus 

H  =  lo^a-e-  +  10^tr»  ^  +  lO^cS-  ^  +  10^a»  ^  -  lC?e  +  H', 
”5“  -y-  ~y-  (^2) 


Where  H'  Is  the  constant  of  Integration. 

Dividing  Equation  (A2)  by  T  =  10^  gives  the  enthalpy 
per  degree.  It  Is  convenient  to  define 


H-Ho  =  a-»-t>©-  +c^^  +  d^^-e'0-"^+H'  -  Ho 

T-  -r  “zr 


(A3) 


Where  the  subscript  o  refers  to  the  thermodynamic  state  at  one 
atmosphere  pressure  and  absolute  zero  of  temperature.  The 
Integration  constant  used  Is  an  average  value  determined  from 


H' 


1  = 


(a4) 


where  Hj^(Tj^)  Is  the  1  “  integration  constant,  obtained  by 

substituting  the  1^^  temperature  and  enthalpy  into  Equation 
(A2)  and  n  Is  the  number  of  such  substitutions. 


For  a  reversible  process  at  constant  pressure  the  entropy 
change  Is 

dS  = 


where  S  Is  the  entropy.  Using  Equation  (Al)  and  again  taking 
the  indefinite  Integral  gives 


Al 
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S  =  alne*  +  be-  +  ^  ^  +  de-  ^  ^  +  S'  ,  (a6) 

2  3  2 

where  the  integration  consstant  S'  Is  determined  in  the  same  way 
as  H'. 


Using  the  Gibbs  free  energy 

G  =  H  -  TS 

it  is  convenient  to  tabulate 

-(0  .  aj=  s  -(h  -  hJ 

T  T 


(A7) 


where  S^  has  been  arbitrarily  chosen  as  zero. 

Using  Equations  (l),(A3)  and  (A5)  for  (K  -  Hq)/T  and  S 

respectively  made  the  machine  computations  unnecessarily  cumber 
some.  Simple  algebraic  manipulations  (nesting)  resulted  in  the 
more  economical  forms. 


C  =  e/-&  ^  b4--&-{c+d-6-) 


H  -  H  =  1 

H'  -  e 

a  -F 

h  +-^  (c/3  -t-  d^/4) 

T  ^ 

105 

2  '  ' 
mm 

S=S'-eO_“  -fa  Ine-  -i-  -e- 
2 


b  +  ’e’{c/2  +  -e-dj 


(A8) 

(A9) 

(AlO) 


A2 


1 
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APPENDIX  B 
Fortran  Listing 

For  those  who  wish  to  use  it,  the  Fortran  source  program 
is  given  here.  It  is  necessary  only  to  punch  the  instructions 
on  Hollerith  cards,  ignoring  the  blanks  Inserted  before  each 
Option. 

Five  comments  occur  in  the  listing.  They  are  marked  by  a 
"C"  in  column  one,  and  are  in  correct  format.  They  will  not 
interfere  with  the  operation  of  the  program,  and  are  inserted 
to  emphasize  the  Options.  On  page  5  the  letter  B  occurring  in 
column  one  identifies  Boolean  instructions.  All  others  are  in 
Fortran II  language. 

No  control  instructions  are  shown  since  they  depend  upon 
the  system  being  used  at  a  particular  Installation. 
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rhtkMODYNAMIC  DAT A  PkObkAM 


CHOOSE  AN  OPTION 

SXGAS  SPD(300»34)  »FORM,(40)»NRtF(300»4)  »CP(5)*A(5»5)»C(5) 

SXLIQ  ALF  L 
21  REAOI.NOPT 

1  FORMAT! 12) 

IF!NOPT)19»30»2 

2  GO  TO  (3»13.i»3f  1000)  iNOPT 

1000  GO  TO  3 

2  IF( SPD( 1 »1 )  )3C01 »3000*3001 

3000  READ  INPUT  TAPE  7 ,4 . ( SPD ( 1 , 1 ) , | FORM ( J ) , J= l  ,40 )  ,SPD ( 1 ,2  ) , 
1SPD(I,3),(NREF(I,<),K  =  1,4),Is1,40) 

A  FORMATtF3.0»40Fl«0»F10.0»2Il*3I4) 

DO  1001  1=41,300 

IS^^n  = 

DO  1001  J=l,40 

1001  SPD  ( 1 ,2  )=SPD(.1 ,2  )+FORM(  J)#SPD(  J»2) 

READ  INPUT  TAPE  7» 1002 ,  ( ( SPD ( I , J ) , J=4 ,34 ) , I  =  1 , 300 ) 

1002  FORMAT! 5F10.0/5F10.0/7F10.0/7F10.0/7F10.0) 

REWIND  7 

3001  GO  TO  ! 6 ,6,22 ) ,NOPT 


C  BEGIN  OPTION  1 

6  READ  7»M0,MDAY,NYEAR 

7  FORMAT! 12,213) 

PRINT  8 ,M0,NDAY,NYEAR 

8  F0RMAT!111H1LISTING  OF  SPECIE-DATa  INPUT  !VAPE7)  TO  THEORETICAL  PE 

IRFORMANCE  PROGRAM  i 

1H/I2)  I2,1H/I2,1 

PRINT  9 

9  FORMAT!103H1  B,p,  .  p 

1  HF!298)  !KCAL/KMOL)  REFERENCES/84H 

^  IKCAL/KMOL)  |G)  (L) 

1  !S)/22H  SPECIE  ID, NO.  MQL.WT.) 

NPAGE=0 

DO  11  1=1,300 

ID=SPD!I,1) 

PRINT  10 , I D,SPD ! I ,2 ) ,SPD! I ,7 ) »SPD ! I ,8 ) , 

1SPD!I,12),SPD!I,13),SPD!I,4),SPD!I»5)»SPD!I,6),!NREF!I,J),J=1,4) 

10  ^FORMAT!9H  1 3  ,  F9.3,2F8. 1  ,SF10.0 ,3H  1 1  I  3 ,2H  1  1 3  !  I 

GALL  SYMBOL!!, 2, NCH) 

NPAGE=NPAGE+1 
IF!NPAGE~50)11, 1005, 1005 
1005  NPAGE=0 
PRINT  9 
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11  CGHTIi'iUE  .  ^ 

PRINT  1007 

1007  F0RMAT(83H1  H(298)-H(0) 

1  C0EFFICIENTS/102H  SPECIE  ID. NO.  PHA.SE  ( KCAL/KROL 

1)  A  B  C  D  E  K1 

1K2) 

M  =  0 

NPAGE=0 

DO  lOOa  I=lt300 
TEST  G=0. 

TEST  L=0. 

TEST  S=0« 

ID=SPD(  III) 

DC  1009  K=li7 
TESTG=SPD( I iK+131+TESTG 
TESTL  =  SPD(  I  iK+ZO+TESTL 

1009  TESTS=SPD( I iK+27)+TESTS 
IF(NPAGE“47 )1021»1020il020 

1020  NPA6E=0 
PRINT  1007 

1021  IF(TESTG)1012il010»1012 

1010  IF(TESTL)1014ilCllfl014 

1011  IF(TESTS) 1016il018f 1016 

1012  PRINT  1013iIOjSPD(I*9)»iSPDlliJ+13)fJ=lf7) 

1013  FORMAT(9H  I3»15H  (G)  F5.0il2H  i 

15F9.5iF8.4iF9. 4) 

NPAGE=NPAGE+1 

CALL  SYMBOL  (Ii2»NCH) 

GO  TO  1010 

1014  PRINT  1015iIDfSPD(i»10)»{SPD(I »J+20) »J=1»7) 

1015  FORMAT(9H  I3»15H  (L)  F5.0»12H  » 

15F9.5iF8.4iF9. 4) 

NPAGE=NPAGE+1 

CALL  SYMBOL  ( 1 iZiNCH) 

60  TO  1011 

1016  PRINT  1017fIDfSPD(Itll)»(SPD(I »J+27)iJ=li7) 

1017  F0RMAT(9H  I3fl5H  (S)  F5.0il2H  i 

15F9,5iF8.4iF9.4i 

NPAGE=NPAGE+1 
CALL  SYMBOL! I f2»NCH) 

GO  TO  1018 

1018  IFCTESTG+TESTL+TESTS)  1019 il008  il019 

1019  N=N+1 

1008  CONTINUE 
PRINT  2OOO1N 

2000  FORMAT (22H0DATA  IS  AVAILABLE  FOR  l4i9H  SPECIES.) 

GO  TO  31 


C  BEGIN  OPTION  2 

13  read  14 

14  FORMAT! 37H  COEFFICIENT  DETERMINATION  FOR  ) 

READ  17»lCP(I)»Isl«5)»H,S 

17  FORMAT(7F10eO) 
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A(1,1)=-3«2040805 
A(  1»2)  =  10. 955440 
A(l,3)=-12. 815569 
A(1,4)=8-.7032594 
Ad, 5)=-2. 6390499 
A(25l)=l. 9376145 
A(2,2l=-7. 2087899 
A(2,3)=10. 775001 
A(2,4)=-8. 0360433 
A(2,5  )=2, 5322178 
A(3,l)=-0i41001527 
A(3,2)=l. 6033583 
A(3,31 =-2i7948590 
A(3,4)=2. 4161851 
A{3,5)=-0. 81466908 
A(4,l)=0. 029156571 
A(4, 2)=-, 11770102 
A(4, 31=0.22712088 
A(4,4)=-0. 22307338 
A(4, 51=0.08449696 
A(5, 11=2. 6473247 
A(5,2l=-5. 2323073 
A(5, 31=4. 6083062 

A('5,4l=-2. 8603279 
A(5, 51=0. 8370043 
DO  15  I=lf5 


15 


16 

117 

18 


C(I1=0. 

DO  15  J=l»5 

C( I )=A{ I ,J)*CP(J1+C( I ) 

XK1  =  3. «-H-3. *C(  1 1-4. 5«C( 2  )-9.*C(3  1-20. 25*0(41+0(51/3. 
XK2=S-1.0986 1*0(1 1-3. *0(2 1-4.5*0(31-9. *0(4 1+0(5 1/18. 
PRINT  16 
FORMAT (IHl) 

PRINT  14 


print  18,(CP(I1 »I=1,51,H,S»(C(I)»I=1»5)»XK1»XK2 


FORMAT! 7 IHO  MATCH  TEMPERATURES 

1  4500  5700/27H 

10C0)-H{01 1/3000  =F8.3/25H 
1  A  B 

IE  '  :<1  K2/12H 

GO  TO  31 


1000  1600  3000 

CP  5F9.3/28HO  (H(3 

5(30001  =F8.3/105H0 
0  D 

1P7E14.51 


C  BEGIN  OPTION  -1 

19  READ  20»TZ»TF»DT 

20  FORMAT! 3F10.0) 

PRINT  21 

21  format (45H1THERMODYNAMIC  FUNCTIONS  FOR  THE  ABOVE  SPECIE) 
GO  TO  26 
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c 

22 

23 


B 

3 


1022 

1023 

1025 

1024 

1026 

24 


1027 
1030 

1029 

1028 


1031 

1032 

1033 

1034 


26 

27 


28 

29 

30 


BEGIN  OPTION  3 
READ  23»NSP»PHASE»TZ*TFiDT 
FORMAT! I3»1A1»3F10.0) 

TESTG= ( PHASE+XGAS )* ( -I PHASE^XGAS ) ) 

TESTL= ( PHASE+XLIQ )* (-! PHASE*XL IQ ) ) 

IF(TESTG)  1022»1023*1022 
IF(TESTL1  1024»1025»1024 
J  =  13 

GO  TO  1026 
J  =  20 

GO  TO  1026 
J  =  27 

DO  24  I=li5 
K=I+J 

C( I )=SPDINSP»K) 

XKl=SPD(NSP»J+6) 

XK2=SPD(NSP»J+7) 

IF(TESTG)  1027»1030.1027 

IF(TESTL)  1028*1029,1028 

PRINT  1031 

GO  TO  1034 

PRINT  1032 

GO  TO  1034 

PRINT  1033 

GO  TO  1034 

FORMAT! 39H1THERMODYNAMIC  FUNCTIONS  FOR  IG)) 

F0RMAT!39H1THERM0DYNAMIC  FUNCTIONS  FOR  !L)) 

FORMAT! 39H1THERM0DYNAMIC  FUNCTIONS  FOR  IS)) 

CALL  SYMBOL!NSP*30»NCH) 


GO  TO  26 

PRINT  27 

FORMAT I97H0 

T 

CP 

IH-HIO) : 

l/T 

S 

IG-GIO) 

)/T)/97H 

!K) 

IKCAL/KMOL-K) 

IKCAL/KMOL-K) 

IXCAL/KMOL-K) 

!KC 

1AL/KM0L-K)/1H  ) 

NSTO?=l.+(TF-TZ)/DT 
DO  28  I=l*NSTOP 
XI  =  I 

T=TZ+!XI-1. )*DT 
THETA=T/1000. 

CPZ=C!l)+THETA#!C!2)  +  THETA«!C!3)+THETA-itC!4)  )  ) +G !  5  ) /TH,ETA**2 
HOT=Xi<l/THETA+C!l  )  +  .5*C!2)*THETA+.33333333*C!3)*THETA**2  +  .25 
1*C!  4)*THETA^f*3“C!5) /THETA**2 

SZ=XX2+C! l)*LOGEF! THETA )+C!2)*TriETA+,5*CI3)*THETA**2+.33333333 
1*C!4)*THETA**3-C!5) /2,*THETA**2 
GZ=HOT-SZ 

PRINT  29*T*CPZfHOT,SZ,GZ 
FORMATI6H  F6.0*3H  1P4E20.5) 

GO  TO  31 
CALL  ENDJOB 
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